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atomic Fermi gas. We develop estimates of spin and thermal diffusivities and discuss magnetocaloric 
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signatures of the spin Seebeck effect. We study an exactly solvable model that illustrates the role 
of momentum-dependent scattering in the magnetocaloric effects. 
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I. INTRODUCTION 

The transport properties of condensed matter systems 
are often measured by driving currents externally and 
measuring the resulting voltages or temperature differ- 
ences. In cold atomic gas clouds, on the other hand, 
transport is more often measured by setting up transient 
out-of-equilibrium initial conditions and measuring the 
subsequent relaxation towards equilibrium [1-5]. In the 
approximation that the cloud is isolated and has an in- 
finite lifetime (no loss of atoms or exchange of energy 
with any degrees of freedom outside of the gas cloud), 
the conserved currents of interest include the energy cur- 
rent and currents of each of the atomic species present. 
In the absence of optical lattices or random potentials 
that violate momentum conservation, one can also ask 
about the transport of momentum (viscosity). 

In this paper we consider the diffusive transport of 
heat and of atoms. We mostly focus on the case of a 
two-species Fermi gas with only inter-species contact in- 
teractions, but start with a somewhat more general dis- 
cussion here. The system may, in addition to diffusive 
transport, also have underdamped or propagating sound 
or other "collective" modes. A gas cloud in a smooth trap 
will have such sound modes, with the longest- wavelength 
sound modes being the often-discussed collective modes 
of the cloud's oscillations within the trap. Here we con- 
sider a gas cloud in a smooth trap, with the cloud at 
global mechanical equilibrium, so that any pressure gra- 
dients in the cloud are sufficiently balanced by trapping 
forces that no underdamped sound or collective modes 
are excited. We also assume that the cloud is everywhere 
near local thermodynamic equilibrium, so the local tem- 
perature T(r) and local chemical potentials Hi{v) can be 
defined. However, the cloud may still have gradients in 
the local temperature and in the local chemical potentials 
of the various species of atoms. If the equilibrium equa- 
tion of state of the system is known (for the unitary Fermi 
gas, see [6, 7]) then measurements of the local densities 
of each species allows these gradients of T and the /ij's to 
be measured. Thus, for example, the local densities can 
be used as local thermometers to allow a measurement 
of the thermal diffusivity by an approach similar to that 



used in [1, 2] to measure the spin diffusivity (but with 
an initial temperature gradient instead of a composition 
gradient). 

The transport currents that we examine in this paper 
are those that arise in linear response to these gradients. 
In a trap, convection currents may also appear in linear 
response, as temperature and/or composition gradients 
may produce density inhomogeneities, and the "heavier" 
regions of the gas cloud will sink towards the bottom 
of the trap while the "lighter" regions rise. These con- 
vection currents are damped by the viscosity. Convec- 
tion will be strongest in wide clouds and should be much 
weaker in high-aspect-ratio clouds with the gradients in 
T and the /ij's oriented along the long axis of the cloud. 
In most of this paper, for simplicity we consider a gas 
in a spatially uniform potential, so such convection cur- 
rents do not appear in linear response. Then mechanical 
equilibrium is indeed a sufficient condition to have a con- 
vectionless gas [8]. 

Quite generally, a temperature gradient drives a heat 
current and a gradient of chemical potential difference 
drives a composition ("spin") current, consisting of op- 
posing currents of the two (or more) atomic species. We 
call these "direct" responses to a temperature gradient 
and a gradient of chemical potential difference "primary 
currents" . In addition to these "primary currents" , there 
are the magnetocaloric currents, namely spin Seebeck 
currents (spin currents induced by a temperature gradi- 
ent) and spin Peltier currents (heat currents induced by 
gradients of chemical potential difference). These mag- 
netocaloric effects have been one of the central research 
topics in the field of spintronics [9] . The spin Seebeck ef- 
fect [10, 11] and the spin Peltier effect [12] have already 
been observed in condensed matter systems, while they 
are yet to be detected in cold atomic clouds. In this pa- 
per we discuss the origin and the physics of these effects 
in a cold atomic Fermi gas, and estimate how large these 
effects can be in realistic experiments. Note that Ref. 
[13] discusses a rather different situation that they are 
also calling the "spin Seebeck effect": they consider an 
unpolarized gas with the two species at different temper- 
atures (thus not in local equilibrium) and a gradient in 
this temperature difference. 



2 



This paper begins with a general discussion of the 
two-species universal Fermi gas and an introduction of 
transport coefficients, spin and thermal diffusivities, spin 
Seebeck effect, and spin Peltier effect. Then we present 
rough estimates of diffusivities and determine the signs 
of the spin Seebeck and spin Peltier effects based on 
physical arguments. In the following section, these qual- 
itative descriptions are justified by approximate solu- 
tions of the Boltzmann transport equation in the classical 
regime. We then compute experimentally verifiable sig- 
nals of the spin Seebeck effect. These results are tested 
against an exactly solvable model, namely atoms with a 
"Maxwellian" scattering cross section. We also use the 
Kubo formula to derive various general relations among 
the transport currents and coefficients. 

The various effects discussed in this paper are probably 
most accessible experimentally for the unitary Fermi gas 
at temperatures of order the Fermi temperature, where 
the diffusivities are at their smallest, so the diffusive re- 
laxation towards equilibrium is slowest and most easily 
studied. 

We set Boltzmann's constant fcg = 1 but explicitly 
keep Planck's constant H. 



II. UNIVERSAL FERMI GAS 

The universal Fermi gas is a two-species Fermi gas with 
only contact (s-wave) inter-species interactions that is 
realized to a very good approximation in recent experi- 
ments with ultracold atoms [14, 15]. We consider such 
a gas in three-dimensional space, with the two species 
having equal mass m. There is no optical lattice, only a 
possible smooth trap potential. The interaction is spec- 
ified by the scattering length a, which can be set to any 
value in experiments by tuning through a Feshbach res- 
onance [16, 17]. For a = this is the standard textbook 
nonintcracting Fermi gas, while for weakly attractive a 
it is very close to the model used by Bardeen, Cooper 
and Schrieffer (BCS) to explain superconductivity. In- 
deed, this system shows paired-fermion superfluidity at 
low temperatures. The limit of infinite |a| is the strongly- 
interacting unitary Fermi gas. 

As is standard, we call the majority species with num- 
ber density "up", and the minority species "down", 
n± < n^. The total number density n = n^ + n^, together 
with m and H set the characteristic length, time and en- 
ergy scales. The scaled dimensionless properties of this 
universal Fermi gas then depend on only three dimen- 
sionless parameters, which can be chosen to be \/{kpa), 
T/Tp and the polarization p = (n-j- — n^)/n. We use a 
convention that the Fermi wavenumber and temperature 
kp and Tp are defined by the total density, so that at high 
polarization Tpj, <C Tp s» Tp-f and kp\, <C kp rj kp^. 
This universal Fermi gas has a variety of regimes of be- 
havior : The polarization p can be low or zero so = 
or it can be near one so ri\ <C ru. The temperature can 
be higher, T > Tpf, or lower, T < Tp±, than both Fermi 



temperatures or, for p > it can be in between them, 
Tpi < T < Tp^. The scattering can be near unitarity so 
\kpa\ is of order one or more, or it can be far from unitar- 
ity so \kpa\ <C 1. At high T it also matters whether \a\ is 
larger or smaller than the thermal de Broglie wavelength 
A = ^/2irh 2 /mT ~ T -1 / 2 . There are also important dif- 
ferences between the a < (BCS side of the Feshbach 
resonance) and a > (BEC side) regimes. 

III. TRANSPORT 

The conservation laws of this Fermi gas are: total en- 
ergy (E) , total momentum (II) , and the total number of 
each of the species (N^ and iVj.). The viscosity measures 
the transport of momentum, which we mostly do not con- 
sider here. Thus we consider primarily the transport of 
atoms and of heat. If there is a nonuniform pressure in 
the system that is not balanced by a trapping potential, 
the gas will accelerate and this will produce free motion 
or propagating sound waves. Here we consider the dif- 
fusive spin and heat transport in a gas with no trapping 
potential and spatially uniform pressure P, so it is at 
mechanical equilibrium. The gas is near local thermo- 
dynamic equilibrium, but with possible weak gradients 
in the local temperature and/or the spin polarization. A 
smooth trap potential may be added via the local density 
approximation (LDA). 

In general, an inhomogeneity of the Fermi gas consists 
of gradients in the local temperature and of the local den- 
sities of the two atomic species. Mechanical equilibrium 
imposes a constraint on these gradients and thus there 
are only two independent linear combinations of the three 
gradients. One way of describing the diffusive dynamics 
is in terms of the atomic densities n, and the currents 
ji of each species i =t,l, leaving the temperature and 
the heat current implicit, since they are dictated by the 
equilibrium equation of state, e.g., T(P, rif, nj.). This de- 
scription of the transport has the virtue that it is in terms 
of what appears to be the most accessible local observ- 
ables in experiment, namely the local densities of each 
species. Since the system in the absence of a trapping 
potential is Galilean-invariant, we have a certain amount 
of flexibility in what inertial frame we use to specify the 
currents. For most of this work, we consider the frame 
where the center of mass of the whole cloud is at rest 
and let the gas have long-wavelength temperature, den- 
sity and/or composition modulations, but always with 
a spatially uniform pressure. The diffusive currents are 
related to the density gradients as 

The diffusion matrix in (1) has two eigenmodes. At 
zero polarization, the symmetry between f and J, dictates 
that the eigenmodes are the spin, which is odd under ex- 
changing species, and the even mode is both the energy 
and the total density, since varying T at fixed P means a 



3 



variation also in n, due to the thermal expansion of the 
gas. Away from zero polarization, when ^ , the dif- 
fusive eigenmodes are in general no longer purely spin or 
purely not spin, instead they are mixtures, thus produc- 
ing the spin Seebeck and Peltier effects. The eigenmode 
where the currents of the two species are in opposite di- 
rections we will call the "spin" mode with diffusivity 
while the other mode where they are parallel we will call 
the "thermal" (or heat) mode with diffusivity Dt- 

Another standard representation of the transport ma- 
trix in terms of the heat current jheat arL d spin current 
j spin is the following: 



is the most directly accessible version of the spin See- 
beck coefficient: if one can set up an initial condition at 
mechanical equilibrium and local thermodynamic equi- 
librium with a temperature gradient but no polarization 
gradient and then measure the resulting spin current, this 
is a measurement of Se' and does not require knowledge 
of the equation of state. 

Note that the three representations, Eq. (1), Eq. (2) 
and Eq. (5) are related by the equation of state, the 
mechanical equilibrium condition, and definitions of spin 
and heat currents. Hence, they are equivalent. 



iheat 
ispin 



K Pe 
Se a s 



VT 



(2) 



where n is the thermal conductivity and a s is the spin 
conductivity. Se and Pe are the spin Seebeck and Peltier 
coefficients, respectively, and they are related by the 
Onsager relation, Pe — TSe. This matrix explicitly 
shows the direct responses (diagonal elements) and mag- 
netocaloric effects (off-diagonal elements), and is in the 
form that is given by the Kubo formula, as we discuss 
below. 

The currents in (2) must be defined properly so that 
they are the transport currents, namely the currents of 
heat and spin relative to the average local motion of the 
gas. Let the local current density of atoms be j„ = Jt+Jt- 
These atoms carry the average heat, sT, where s is the 
average entropy per particle, and the average spin polar- 
ization p. Therefore the local heat and spin transport 
currents arc 



jheat = je — MtJt ~~ M-Jj- — S ^jr. 
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where j e is the local energy current. Since the above cur- 
rents measure only the transport relative to the average 
motion of the gas, they are reference frame independent. 
For more details of definitions of currents, see e.g. Ref 
[18]. 

If the full equation of state of the system is known, 
then measurements of the pressure and the local densi- 
ties can be converted to local temperatures and chemical 
potentials. But the local density n and polarization p 
are directly observable without requiring knowledge of 
the equation of state, so yet another convenient form of 
the transport equations is 
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At p — 0, we have k' = K and D s — 2) s , but when 
these quantities in general differ due to the mixing be- 
tween spin and heat transport. It is possible that Se' 



IV. DIFFUSIVITIES 

Let's first present rough "power-counting" estimates 
of the spin and thermal diffusivities, D s and Dt, respec- 
tively. At the level of power-counting the differences be- 
tween the various possible definitions of these diffusivities 
are small and are ignored here. Previous work [1, 19] on 
the unpolarized gas (p — 0) shows that D s for T > Tp 

is the larger of ^(^-) 3 ^ 2 and „_J Sr. In the re- 
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cent experiment, which was performed at unitarity [1], 
this high-T behavior is observed, with significant devi- 
ations apparently beginning between T = 2Tp and Tp 
as Tp is approached from above. Staying in this high-T 
regime, as we move to high polarization (p near 1) at a 
given n and T, the scattering time of species i is roughly 
Tj ~ n 1 <yv where i ^ j and v r = |v^ — vj and a is the 

s— wave scattering cross section (see Eq. (12)) evaluated 
at a typical value of momentum. Thus, the scattering 
time of the down atoms decreases by only a factor of 
two due to the increase of the density n+ of the up atoms 
that they scatter from. The up atom scattering time r-j-, 
on the other hand, increases by a factor of /nj, as the 
down atoms that they scatter from become dilute. At 
high polarization, the spin current consists of the down 
atoms moving with respect to the up atoms at typical 
speed v± ~ ^/T/m, so D s ~ wfrj, is not strongly polar- 
ization dependent for T > Tp^; the experimental results 
[1, 2] are consistent with this. The heat, on the other 
hand, is mostly carried by the up atoms at high polar- 
ization, resulting in Dt ~ D s n^/n±, a relation between 
the two diffusivities that appears to remain true at high 
polarization for all T away from the superfluid phases. 
At p — and high T the two diffusivities are compa- 
rable, but Dt remains larger than D s because a single 
s-wave scattering event completely randomizes the total 
spin current carried by the two atoms, while the compo- 
nent of the heat current carried by their center of mass 
is preserved. Thus it appears that the heat mode always 
diffuses faster than the spin mode. 

Moving towards lower T, let's next pause at T — Tp^, 
noting that here ~ ~ y/T/m, and is the larger 
of ^ and Tk 2 a 2 ■ For the polarized gas ~ r^n^/n^. 

We next (still just power-counting) look at the polar- 
ized gas in the intermediate temperature regime Tpt < 
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T < Tp^ where the majority atoms are degenerate 
(v-f y/Tp^Jm), while the minority atoms are not 
(ui ~ yjT/m). Some changes from the high-T regime 
are: only up atoms with energy within ~ T of Tp^ are 
involved in the scattering and all but a fraction T /Tp^ 
of the final states of the scattering are Pauli-blocked due 
to the degeneracy of the up atoms. This increases 
by a factor of (Tp^/T) 2 , and r-f by a factor of T F ^/T, 
compared to their values at T = Tp^. Thus, is the 

larger of ^ (3p) and Tp ^ a2 (^f 1 ) . As a result, 

D s ~ vfr^ is the larger of and m ^P a 2 , while D T is 

again larger than D s by a factor n^/n^. This estimate of 
-D s is consistent with a previous quantitative calculation 
[20, 21]. Note that the temperature dependence of D s 
crosses over from a decreasing function of T at low T to 
an increasing function at high T. The recent measure- 
ments [2] of the spin drag in a polarized unitary gas show 
this crossover occurring at roughly T = QATp^. On the 
BEC side of the Feshbach resonance, the minority atoms 
bind in to bosonic Feshbach molecules at low enough T. 
But as long as these molecules remain nondegenerate and 
thus not superfluid, the above estimates of the diffusivi- 
ties should hold. 



For T < Tpj,, the minority atoms become degenerate. 
This leads to superfluidity on the BEC side of the Fesh- 
bach resonance as well as at low polarization near unitar- 
ity. But there are regimes on the BCS side of the reso- 
nance as well as near unitarity at high polarization where 
the minority atoms (near unitarity strongly "dressed" as 
polarons) form a degenerate Fermi gas. Here the impor- 
tant change from the intermediate temperature regime 
at the level of "power-counting" is that only minority 
atoms with energy within ~ T of Tp\. are involved in 
the scattering and they have momentum kp\ instead of 
a thermal momentum. This increases the diffusivities by 
a factor of Tp^/T, so the spin diffusivity in these degen- 
erate Fermi liquid regimes is the larger of ^"^ff^ and 
mT^l^a 2 ' ^ e ex P ec t that Dp is still greater than D s by 
a factor of n^/n^ but this question should be examined 
more carefully within Fermi liquid theory. 

At low temperature, the polarized Fermi liquid may be- 
come a p-wave superfluid, with pairing within one species 
mediated by the attraction to the other species [22]. Or 
it may have a FFLO phase with Cooper pairs of nonzero 
total momentum [23-25]. In the superfluid phases, the 
thermal diffusivity should diverge to infinity, as heat is 
carried ballistically by second sound modes. The spin 
diffusivity presumably remains finite in the superfluid 
phases, although, as we discuss below, the spin Seebeck 
effect appears to generally be divergent in a polarized 
superfluid. 



V. SPIN SEEBECK AND SPIN PELTIER 
EFFECTS 

More challenging to estimate than these spin and ther- 
mal diffusivities are the effects that mix spin and heat 
transport, namely the spin Seebeck and the spin Peltier 
effects. Here we present "simple" arguments for the signs 
of these effects in two regimes: (1) far away from uni- 
tarity for all temperature ranges, and (2) at unitarity 
in the classical regime. Here we always consider a spin- 
polarized gas, since these "magnetocaloric" effects vanish 
by symmetry in the case of an unpolarized gas where the 
two species also have equal mass. 

(1) Well away from unitarity (|&f&| "C 
mm{l, y/Tp/T}), the scattering cross section is es- 
sentially a 2 and independent of momentum for all 
temperatures, since the interaction is weak. Generally, 
the scattering rate is proportional to the cross section 
x relative speed. Thus, in this regime, the scattering 
rate is ~ a 2 |v^ — vjj. This implies minority atoms 
will scatter more frequently with the majority atoms of 
higher energy than majority atoms of lower energy since 
higher-energy majority atoms have higher relative speed. 
Consequently, the direction of the minority current is 
along the flow of higher-energy ("hot") majority atoms 
and thus parallel to the heat current. For a uniformly 
spin-polarized gas with a temperature gradient, where 
the spin Seebeck effect occurs, the primary current is the 
heat current transporting "hot" majority atoms from 
the hot region to the cold region and transporting "cold" 
majority atoms from the cold region to the hot region. 
The fact that the minority current is aligned with that 
of the "hot" majority atoms means the direction of 
the net spin current (spin Seebeck current) is opposite 
from that of the heat current. Thus, the initially cold 
region becomes less polarized due to minority atoms 
transported by the spin Seebeck current. 

For a polarized gas with a polarization gradient but 
zero temperature gradient, where the spin Peltier effect 
occurs, the primary current is the spin current which 
transports minority atoms from the less-polarized region 
to the more-polarized region. Since again these minor- 
ity atoms scatter more often with "hot" majority atoms, 
the resulting heat current is towards the more-polarized 
region, resulting in a spin Peltier (heat) current whose 
direction is opposite to the primary spin current. In sum- 
mary, for a gas far from unitarity, the primary currents 
and the magnetocaloric currents are in opposite direc- 
tions. In other words, the off-diagonal elements in Eq. 
(5) are negative while the diagonal elements are positive. 

(2) At high temperatures (T 3> Tp^) and near uni- 
tarity (|fejra| 3> kp\ > 1), the s-wave scattering cross 
section is l/(k r /2) 2 , where k r — |k^ — kj,| is the rela- 
tive momentum. Therefore, the scattering rate is roughly 
~ ^2 x v r ~ . As a result, now minority atoms scat- 
ter more often with "cold" majority atoms. Since this is 
exactly the opposite from the case of far away from uni- 
tarity, the spin Seebeck and spin Peltier currents are re- 
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(a) Away from Unitarity: 

Scattering rate ~ <r|v r — VjJ ~ a 2 |vj — Vj | 



(b) Near Unitarity: 

Scattering rate ~ it|vj — vj 
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FIG. 1: (Color online) Illustration of the spin Seebeck effect, 
(a) Far away from unitarity, the scattering rate is proportional 
to the relative speed. Thus, "hot" majority (f) atoms (top) 
collide more often with minority (J,) atoms than do "cold" 
majority atoms (bottom), giving the spin Seebeck current and 
the heat current opposite directions, (b) At unitarity, the 
scattering rate is inversely proportional to the relative speed. 
Therefore, the direction of the spin Seebeck current is reversed 
relative to (a). 



forward to apply the above argument to predict the di- 
rection of spin Seebeck and/or spin Peltier currents since 
the many body effects may significantly modify the scat- 
tering cross section [26, 27], which begins to depend on 
the center of mass momentum as well as the relative mo- 
mentum. There is, however, a different line of argument 
that indicates that the sign of the spin Seebeck effect near 
unitarity remains the same as the temperature is lowered. 
Consider low enough temperatures and polarization less 
than the Chandrasekhar-Clogston limit [28, 29], in the 
superfluid phase [30, 31]. When there is a temperature 
gradient in the system, heat flows "ballistically" from the 
hot region to the cold region by flow of the normal fluid 
with respect to the superfluid (in the usual two-fluid de- 
scription of the superfluid phase) . In the reference frame 
where the total particle density current vanishes (center 
of mass frame), the superfluid flows in the opposite di- 
rection to counterbalance the mass current of the normal 
fluid. Since the s-wave superfluid consists of equal num- 
bers of majority and minority atoms, both the spin cur- 
rent and the heat current are carried only by the normal 
fluid. As a result, the spin Seebeck current and the heat 
current are in the same directions at low temperature in 
and, presumably, near the superfluid phase. Therefore, 
we expect the sign of the spin Seebeck effect to remain the 
same for all temperatures at unitarity. Both the thermal 
conductivity and the spin Seebeck coefficient will diverge 
at the transition to the superfluid phase. 

As another approach to these questions, there is an 
interesting artificial interaction, namely the Maxwellian 
interaction where the scattering cross section is pro- 
portional to l/v r and the Boltzmann equation can be 
solved exactly in the high temperature limit. Since 
the scattering rate does not depend on relative veloc- 
ity ((l/u r ) x v r — constant), there are no spin Seebeck 
or spin Peltier currents in this case. Then, if we perturb 
the scattering cross section around the Maxwellian case, 
putting in an additional relative-velocity dependence to 
the scattering rate "by hand" , we can perturbatively cal- 
culate the spin Seebeck and spin Peltier currents and 
manipulate the direction of these currents by changing 
the sign of the perturbation. This allows us to explicitly 
show how the magnetocaloric currents are generated from 
a relative-velocity-dependent cross section. This will be 
discussed in the following section in detail. 



versed relative to the above discussion in (1). Therefore, 
at high temperature and unitarity, the primary currents 
and the magnetocaloric currents are in the same direc- 
tions, giving Se' in Eq. (5) positive sign. As we will see 
in the next section, the spin Peltier coefficient Pe' in Eq. 
(5) is negative for low polarization and becomes positive 
for high polarization. This sign change in Pe' comes from 
the definition of currents and choice of driving forces and 
this will be clarified in the section VII where we discuss 
the Kubo approach. Directions of the spin Seebeck effect 
in both limiting regimes are illustrated in figure 1. 
At low temperatures and unitarity, it is not straight - 



VI. QUANTITATIVE APPROACH 1 - THE 
BOLTZMANN EQUATION 

A. Linearized Boltzmann equation and its scaling 

In the limit of high temperature T ^> Tp, the gas 
is effectively classical, and its dynamics obey the Boltz- 
mann equation. In the absence of external forces but 
with gradients in local temperature and local densities, 
the steady-state Boltzmann equation for species i (i =j" 
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or I) is the following (i ^ j) 

h. 



V/, 
d 3 k 



%da 



fi|k t 



^(/(k')/(k')-/(k t )/(k ; )), 



(27T) 

(6) 

where /(kj) is the momentum distribution of species i, 
and the velocity is Vj = ftk,/m. kj. and k^ are mo- 
menta after collision and satisfy the energy momentum 
conservation. Working near equilibrium, we linearize the 
Boltzmann equation by introducing a small deviation fy: 

/i = /?(l + Vi), (7) 

where is the equilibrium distribution. In this high 
temperature regime, the equilibrium distribution is the 
Boltzmann distribution, ff{k) — t^A 3 exp(— E{k)/T) 

and E(k) = *t£t- Since the most relevant length scale in 
the classical regime is the thermal de Broglie wavelength 
A, it is convenient to scale the wave vector k with A; k 
= q/A. Then, E{q)/T = £. 

Let's impose the mechanical equilibrium condition. In 
this classical limit, it is enough to use the ideal gas pres- 
sure at equilibrium; P = nT = (n^ + n{)T. From the 
spatially uniform pressure condition we can relate the 
density gradients and the temperature gradient via 

to linear order in the gradients. This relation enables 
us to express the currents in terms of any two linearly- 
independent "driving" terms such as (Vn^,Vn^), 
(VT, riVp), (VT, V(/i| — hi)) or any other convenient 
combinations. It is convenient to work with (Vn-j-, Vn|) 
in intermediate stages of the calculation and then trans- 
form it to the desired combination of driving forces using 
the equilibrium equation of state. Then the A scaled lin- 
earized Boltzmann equations we need to solve become 
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with i ^ j. ifi'i takes q^ as an argument. 

We follow the standard definitions of the particle and 
energy currents of each species [32] , 
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(11) 

Because of Galilean invariance, we need to specify an 
inertial frame. Except when specified otherwise, we work 
in the frame where the local particle current is zero: j n 
0. 



B. Approximate Solution 

The s-wave scattering cross section that captures most 
physics of a short-range interaction is 
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where we scaled the scattering length a with A. An ex- 
act solution of the Boltzmann equation for such a cross 
section is not known and thus we need to resort to ap- 
proximation methods. One of the standard ways to find 
an approximate solution of a linearized Boltzmann equa- 
tion is the moment expansion method (for example, see 
[33]). Considering all symmetries and assuming the true 
solution is analytic near small q, we take the following 
ansatz for ipi (i =f,J,): 
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We only consider the case where Vri| and Vn^ are both 
parallel to the z axis. We need to determine the dimen- 
sionless coefficients {Cuj}- 

The procedure to obtain an approximate solution of 
the Boltzmann equation is the following: First, insert the 
above ansatz into the right hand side of Eq. (9). Then 
multiply both sides of Eq. (9) by qf e qi Z j^js ex P[ — if /47r] 
(I = 0, 1, 2, ...L) and integrate out all momenta. Match- 
ing coefficients of density gradients gives 4(L + 1) lin- 
ear equations for the {Cuj}, all of which, however, are 
not linearly independent due to Galilean invariance. We 
need to fix the reference frame to uniquely determine a 
solution. Once we choose an appropriate reference frame 
(usually j-j- + ii = 0), we have 4(L + 1) linearly indepen- 
dent equations for the {Cuj}. Determining the {C^j}, 
we have an approximate solution to the Boltzmann equa- 
tion, i.e. an approximate momentum distribution from 
which we can calculate all currents of interest. Here we 
present results of two limiting cases, far away from uni- 
tarity (j^ — a 2 ) and at unitarity (J^ = p.), which al- 
low an analytic solution (of this approximation) without 
special functions. These correspond to the two limits 
A/|a| ^> 1 and \/\a\ <C 1, respectively. For a general 
scattering length a, it is still possible to find a closed 
form expression in terms of exponential integrals and in- 
complete Gamma functions whose arguments depend on 
A/14 

Obtaining Cuj from a straightforward calculation in 
the j-j- + jj, = frame, we can express the heat current 
and spin current in terms of Vn^ and Vn^. Here we 
choose to express final results in the format of Eq. (5) 
since we want to study the spin Seebeck coefficient Se' in 
detail, which could be the most directly accessible signa- 
ture of the spin Seebeck effect in experiments. Therefore, 
we transform these two gradients to VT and nVp, us- 
ing the equilibrium equation of state and the mechanical 
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equilibrium condition. The transformation matrix is 

Since the left hand side of Eq. (9) is a third-order poly- 
nomial in q, the simplest ansatz is with L — 1. In princi- 
ple, we can go up to any order in L we want, but an L > 2 



ansatz complicates the computation, while the simplest 
ansatz already exhibits nontrivial results. Furthermore, 
we find that the change in Se' on moving from the L = 1 
to the L = 2 approximation is quite small: about a 1% 
change near unitarity and in the value of A/|a| at the zero 
crossing, growing to near 7% far from unitarity. There- 
fore, here we only present results of L = 1. 



We present our results in the conventional format (without A scaling): Near unitarity (\a\ ^> A), 

J h eat\_ 45**/' (h ( T y/ 2 \ / < \ ( ^^-f Tln(^)\ /VT\ 
j sp in/ 608\/2 V^W y \ ( "t-"4.) |2 J [nVpJ 



5 / (73n2 + %2r^nx + 73n?) 
n 1 



16 V n-j-nj. 
Far away from unitarity (|a| <C A), 

jspin/ 224^2 ym(k F a) 2 \l T F J y ("t~"+) 43 y ^nV^y 

where k' q (proportional to the thermal conductivity far away from unitarity) is 

. 5 A29n? + 26nt% + 29n?)\ / % 



(15) 



where (proportional to the thermal conductivity at unitarity) is 



(17) 



+ (n T -n ; )ln I -± . (18) 



The above matrices clearly exhibit the existence of the 
spin Sccbeck and spin Peltier effects (non-vanishing off 
diagonal terms) only for nonzero polarization. We will 
mostly focus now on the spin Seebeck coefficient Se', 
which gives the spin current due to a temperature gradi- 
ent in the absence of a spin polarization gradient. 

As argued in the previous section, the spin See- 
beck coefficient changes sign as a function of interaction 
strength. Near unitarity (Eq. 15), it is positive so the 
spin Seebeck current and the heat current are in the same 
direction. Far away from unitarity (Eq. 17), it is nega- 
tive so the spin Seebeck current and the heat current are 
in the opposite direction. 

Next let's consider the polarization dependence of the 
transport coefficients. The heat current and the temper- 
ature gradient are even under spin index exchange while 
the spin current and the polarization gradient are odd. 
Therefore, thermal conductivity and spin diffusivity are 
even functions of polarization while magnetocaloric ef- 
fects are odd functions of polarization. The above matri- 
ces satisfy these polarization parity constraints and the 
form of the polarization dependence of the transport co- 
efficients is the same in both limits of large and small 
A/|<z|. In fact, we can show that the polarization de- 



pendence (thus and dependence) of the transport 
coefficients maintains this form for all values of A/|a| and 
to all orders of approximation. The proof is given in Ap- 
pendix A. Here we study the Seebeck coefficient Se' in 
detail, which is our prime interest. 

Although it is conventional to scale the scattering 
length a with k F (as we did in the above matrices), it 
is easier to see the structure of the Seebeck coefficient 
in terms of A/|<z| in classical regime. Once we obtain an 
approximate solution of the Boltzmann equation with a 
general X/\a\ and L = 1, we can explicitly show that 



Se' = A^l M A/M) 
m nX 6 1 



(19) 



where h\(x) is a dimensionless function that contains 
Ei(x), the exponential integral, and diverges as A/|a| — ¥ 
co (see Eq. (17)). Since it contains no explicit temper- 
ature or polarization dependence, h\(x) is independent 
of polarization and temperature at this order of approx- 
imation (L = 1). Therefore, once we scale the scatter- 
ing length by A and factor out dimensionful parameters 
and polarization, the dependence of Se' on the scatter- 
ing length is determined by hi(x) and the value of A/|a| 
at which Se' crosses zero is solely determined by the 




( a ) 5a r (1, 0) 



Vial 



FIG. 2: Normalized Seebeck coefficient Se' as a function of 
\/\a\ in the classical regime. The normalization is by a fac- 



tor of 



/ h 1 l+47r(q/A) 2 \ 
\mTX 3 (a/A) 2 J 



(n t -ra 4 ) 



This choice of nor- 



608^2 

malization, which is inversely proportional to a typical value 
of scattering cross section with scaled scattering length a/A, 
gives finite values at both limits and gives the Seebeck coeffi- 
cient 1 at unitarity (Eq. (15)). We can see that Se' changes 
sign near \/\a\ ~ 3.62. Since we factored out all explicit tem- 
perature and polarization dependence in the normalization 
choice and scaled the scattering length, this plot remains the 
same for all temperature and polarization ranges at this order 
of approximation (L — 1). 




-0.1 
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(b ) 6m (t, 0) 




Time (s) 



equation h\{x) = 0, which is independent of tempera- 
ture and polarization. In the L = 1 approximation, the 
zero-crossing value is A/|a| ~ 3.62. Figure 2 is a plot of 
the normalized Se' as a function of A/|o|. In the L = 2 
approximation, we find that the scaling function h\{x) 
slightly changes to h,2(x) and the zero-crossing point re- 
mains at A/|a| ~ 3.62. In Appendix A, we show that the 
structure of Eq. (19) (and other transport coefficients in 
a similar manner) remains to all orders of approximation. 
Therefore, we may conclude that in this classical regime 
the Seebeck coefficient is linearly proportional to the po- 
larization p and inversely proportional to TA 3 , once we 
scale the scattering length by A. 

Note that Eqs. (15) and (17) do not explicitly satisfy 
the Onsager relation and the spin Peltier coefficient Pe' 
is still negative for low polarization even near unitarity. 
These come from the definition of diffusive currents and 
choice of representation and will be discussed in detail 
in the next section in terms of the Kubo formula. For 
now, we will focus on the spin Seebeck coefficient Se' 
which appears to be the most promising candidate of the 
magnetocaloric effects to be detected in experiments. 

From the Einstein relation, we obtain the thermal dif- 
fusivity Dt after dividing k' u and n' a hy Cp = 5n/2, 
the heat capacity per volume at fixed pressure and po- 
larization. These results confirm the "power-counting" 
estimates of diffusivities: In case of an unpolarized gas 

/ s 3/2 

(n-f = ni), D s at unitarity is = 1.1^ J , which 
is consistent with previous work [1, 19]. Also, the ther- 



FIG. 3: (Color online) 8rii(t, 0) normalized by initial deviation 
Sriio as a function of time near unitarity. (a) Majority density 
deviation relaxes monotonically for any polarization, (b) Mi- 
nority density deviation shows nonmonotonic relaxation for 
p > due to the spin Seebeck effect. Here we assume 6 Li 
atoms with T/Tf = 4 and the longitudinal length of the trap 
L — 200/im to set the time scale. 



mal diffusivity does satisfy the inequality, Dt > D s at 
zero polarization. Furthermore, for a highly polarized 
gas (n T > ni), D T ~ ^i D s ~ ^ D s, which is also as 
expected from the "power-counting" estimates. 



C. Estimate of the Spin Seebeck Effect 

The spin Seebeck effect seems to be more accessible 
to experiment than the spin Peltier effect, since the mea- 
surement of spin currents has already been done [1,2] and 
seems more straightforward than measuring heat cur- 
rents. Also, the spin mode diffuses slower than the heat 
mode, so the change in spin polarization produced by the 
spin Seebeck effect will relax slowly, enhancing its de- 
tectability. An initially fully equilibrated spin-polarized 
gas could be heated at one end, producing a tempera- 
ture gradient, and then the resulting spin current could 
be measured if it is large enough. 

Let's make quantitative estimates of signatures of the 
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FIG. 4: (Color online) 8rii(t, 0) normalized by initial deviation 
Sriio as a function of time far away from unitarity. (a) Ma- 
jority density deviation relaxes nonmonotonically for p > 
due to the spin Seebeck effect. Inset figure magnifies the non- 
monotonic part of majority density deviations which are very 
weak compared to the minority density deviations at unitarity 
shown in Fig. 3. (b) Minority density deviation relaxes mono- 
tonically for any polarization. Here we assume 6 Li atoms with 
T/Tf = 4 and \/\a\ = 4 and the longitudinal length of the 
trap L = 200/im to set the time scale. 



spin Seebeck effect that are relevant to such a proposed 
experiment. We will make our estimates for a gas in a 
uniform potential, but the results should be roughly cor- 
rect for a gas cloud in a trap if one compares points at 
opposite ends of the cloud that are at the same potential 
so will have the same local densities and polarization at 
equilibrium. First, apply a small, long wavelength tem- 
perature inhomogeneity along the z axis. In mechanical 
equilibrium, nonuniform temperature implies nonuniform 
total density (by Eq. (8)), thus temperature modulation 
implies density modulation. This enables us to write the 
initial total density as 



n(t = 0, z) = no + 5uq cos wz, 



(20) 



where w (= tt/L) is the wavenumber of the modulation, 
L is the length of the system over which the full temper- 
ature difference is applied, and Sno is the small deviation 
of total density from the mean value no due to this tem- 
perature difference (at uniform pressure). Eq. (8) implies 



that if we initially locally heat z = L relative to z = 0, 
then this location initially has lower density because of 
thermal expansion. From the initial condition of uniform 
polarization, the density of each spin component at t = 
is 



rii(t = 0, z) = n i0 + SriiQ cos wz 
n i0 



5n lQ 



7 If) 



Sn . 



(21) 
(22) 



Then, we write a diffusion matrix in the form of Eq.(f), 
assuming that there are only these diffusive currents (no 
initial other motion of the gas). 

As mentioned above in Section III, the diffusion matrix 
has two eigenmodes, the spin mode with the eigenvalue 
D s and the thermal mode with the eigenvalue Dt- Ap- 
plying the continuity equation to the diffusion matrix, we 
obtain a coupled set of heat equations with appropriate 
initial conditions, which can be solved immediately for 
t > in terms of the two eigenmodes. The solution of 
the heat equation gives us the space-time dependence of 
density of each species: 



rif(t,z) 
H{t,z) 



"to 

%0 



-£> 3 u/t 



5n 



cos wz 



(23) 



where a, j3, 7, and £ are determined by initial conditions. 
See Appendix B for details. Since £>t > ®s > 0, the two 
diffusive modes relax to global equilibrium at different 
rates. At nonzero polarization, I7IJCI 7^ 1- From these 
two properties, the density deviations of each species 
have different time evolution from one another: the den- 
sity deviation of one species relaxes nonmonotonically 
while the density of the other species relaxes monoton- 
ically. This nonmonotonic relaxation of density devia- 
tion is a qualitative signature of the spin Seebeck effect. 
Near unitarity, we already know that the heat current 
and the spin current due to a temperature gradient are 
in the same direction. Therefore, the initially colder re- 
gion [z < n/2w) becomes more polarized due to the spin 
Seebeck current. It turns out that the minority density 
deviation changes sign as it relaxes towards equilibrium. 
Far from unitarity, on the other hand, the initially cold 
region becomes less polarized, since the direction of the 
spin current is reversed, resulting in nonmonotonic relax- 
ation of the majority density deviation. Figures 3 and 4 
are plots of relaxation of the density deviations of each 
species as a function of time for the two cases of near 
unitarity and of far away from unitarity. 

Another consequence of the spin Seebeck effect is the 
change in polarization. As argued, the sign of the local 
polarization change depends on the direction of the spin 
Seebeck current. Figure 5 shows the (normalized) devia- 
tion of polarization from the average value as a function 
of time at z = 0. As expected, deviation is positive near 
unitarity and is negative far away from unitarity. 
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FIG. 5: (Color online)Normalized polarization deviation as 
a function of time for three global polarization values, p = 
0.9, 0.6, and 0.3. Near unitarity (solid lines) the polarization 
deviation is positive while far away from unitarity (dashed 
lines) it is negative. We found that the polarization deviation 
is the largest near p = 0.6 in both limits. Same physical 
parameters as Figures 3 and 4 were used to fix the time scale. 
The polarization deviation Sp(t, 0) at the initially cold end of 
the cloud is normalized by the initial total density deviation, 
Sno/rio- 

Perhaps one of the most easily accessible quantities in 
experiment is the density deviation of the species that 
shows nonmonotonic relaxation vs. time. A dimension- 
less measure of the extremum deviation is 3:t ' ' ) , 

Sn i0 

where i =f away from unitarity and i —\. near unitarity 
and tj ex t is the time when density deviation of species 
i is its extremum of opposite sign from the initial con- 
dition. Figure 6 shows the dimensionless measure of the 
extremum density deviation as a function of polarization. 
Near unitarity (solid line), it is a monotonically increas- 
ing function of polarization for p < 1 (at p = 1, nj, is zero 
so the deviation of is undefined). Around p 0.8, the 
strength of the spin Seebeck effect by this measure is 
about 5 %. Far from unitarity (dashed line), this signal 
is much weaker. This is expected since the spin Seebeck 
effect is a consequence of interaction and does not exist 
in the non-interacting gas. Thus, the spin Seebeck effect 
should be strongest near unitarity. Note that in both 
limits, the spin Seebeck effect disappears when the gas is 
unpolarized, p = 0. 

The spin Seebeck effect is a small effect. In the regimes 
where we have been able to estimate it and using the mea- 
sures we have been able to devise, it is less than a 10% 
effect. However, it is worth emphasizing that these com- 
putations are done in the classical regime, so the spin See- 
beck effect does not demand extremely low temperatures 
to detect it. We expect it to be most readily detected 
at temperatures of order Tp, where the diffusivities are 
minimized so the resulting time scales are longest. In 
addition, the procedure to detect it discussed in this sec- 
tion does not require knowledge of the system's equation 



FIG. 6: Extremum density deviation normalized by initial 

deviation, j Sn '^l' n °* t '°' > j as a function of polarization. Near 

unitarity (solid line), i =\,; far from unitarity (dashed line), 
i —f. ti tex t is the time when the density deviation of species 
i reaches its extremum. Far from unitarity (dashed line), the 
normalized extremum density deviation vanishes in the high 
polarization limit since there are no minority atoms to scatter 
from. 

of state. Importantly, we propose a type of experiment 
where the the spin Seebeck effect is a qualitative effect, 
namely a nonmonotonicity of the system's relaxation to 
global equilibrium. 

D. Exactly Solvable Model 

So far, our approach was based on an approximation 
method. Therefore, it is worthwhile to compare the 
main findings to a different approach, namely perturbing 
around the Maxwellian model, where the scattering cross 
section is inversely proportional of the relative speed [8] . 
The linearized Boltzmann equation can be exactly solved 
for the classical Maxwellian model. 

Let's recall the argument from which we determined 
directions of the magnetocaloric effects. Since the scat- 
tering rate is proportional to the product of the cross 
section and relative speed, the scattering rate is inde- 
pendent of momentum for the Maxwellian interaction. 
Therefore, we expect that Se' = for any polarization. 
It is straightforward to exactly solve the Boltzmann equa- 
tion using the ansatz Eq. (13) with L = 1 to confirm this. 

The next step is to perturb the Maxwellian scattering 
cross section to generate magnetocaloric currents. One 
simple way to perturb the cross section is to add a small 
term which depends linearly on relative momentum to 
the original cross section: 

<*> 

where Sq is a constant of the dimension of length and 
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e is a small expansion parameter. For a positive e, the 
momentum dependence of the scattering rate is similar 
to the case away from unitarity, higher scattering rate for 
higher relative momentum. Thus we expect the resulting 
spin Seebeck current is in the opposite direction from the 
primary heat current. For a negative e, the momentum 
dependence resembles the case near unitarity and there- 



fore we expect the resulting spin Seebeck current is in 
the same direction as the primary heat current. 

To find the solution of the Boltzmann equation up to 
linear order in e, we need L — 2 in the ansatz of Eq. 
(13). After a straightforward calculation, we obtain the 
following results: 



f'iheat\ _ _ D ( K M 2n t n^A 2 10 I 1 ^ ^ ) l n ( n\ ) \ ( ^ T \ (25) 

\) spin J ° y ne ^ 2 ~ n ^ _L (1 _ 20p^) J \nVp) ' 

where k! m (proportional to the thermal conductivity of the Maxwellian Model) is 

, 1 / (n| + n T% + n|) 37re (9n| + lOn^ + 9rcJ) \ 7re(nt - fn t \ 



D is m g^ nX 2 , with units of a diffusivity. 

We immediately see that all off-diagonal elements van- 
ish for zero polarization. When e = 0, we see that Se' — 
and thus we conclude that the spin Seebeck effect is a 
consequence of momentum dependence of the scattering 
rate. The sign of the spin Seebeck current at nonzero e is 
as expected. For Pe', we again see the logarithmic term 
which will be discussed in the following section. Perturb- 
ing the Maxwellian scattering cross section re-confirms 
the sign argument for the spin Seebeck effect that we 
presented in the previous Section. 



VII. STRUCTURE OF TRANSPORT 
COEFFICIENTS - KUBO FORMULA 



The Kubo formula gives a formally exact expression 
for the transport coefficients in the linear response regime 
(see e.g. [34, 35]). For irreversible processes in the lin- 
ear response regime, what the the Kubo formula gives 
are the transport coefficients for the dissipative forces 
and currents associated with the entropy production. At 
mechanical equilibrium for our two-species gas, the dis- 
sipative forces are VT and V(/i-f — /i^) [18], thus what 
we obtain from the Kubo formula is the diffusion matrix 
in the form of Eq. (2), whose off-diagonal terms always 
satisfy the Onsager relation. Thus, the diffusion matrix 
in the form of Eq. (5), in which we summarized the re- 
sults in the previous section, does not generally satisfy 
the Onsager relation, although it is possibly easier to ob- 
serve experimentally. In Appendix C, we summarize the 
results in the form of Eq. (2) and explicitly show both 
our approximate solutions and the perturbative solution 
from the exactly solvable Maxwellian model indeed sat- 
isfy the Onsager relation. 

Following Ref. [34] and from Eq. (2), the spin Seebeck 



coefficient Se and Peltier coefficient Pe can be expressed 
via the Kubo formula as 



TSe = 



1 



3VT 
1 



3VT J 



dt(3 heat (0) ■ J 

spin 

(*)> 



(27) 



dt(J heat (t) ■ J 

spin 

(0)) = Pe, (28) 



where V is the total volume of the system and the current 
J is the volume integral of local current density, 



J(t) - J d 3 xj(x,i). 



(29) 



The average is taken over the equilibrium distribution, 
which is just the Boltzmann distribution of each species 
in the high-temperature classical regime. 

Since we assume Galilean invariance, the total momen- 
tum of the entire gas is conserved. Therefore, any physi- 
cal quantity which is transported with the total particle 
current J n remains finite for all time t and gives a di- 
vergent contribution in the Kubo formula, meaning that 
quantity moves "ballistically" rather than diffusively. 
Therefore, it is crucial to use the frame-independent def- 
inition of the local diffusive currents from Eq. (3) and 
Eq. (4). Let's write them again and slightly manipulate 
the heat current: 



_ H 

ispin — Jt 

n 



— J; 

n 



iheat 



5 



- sT 'h 



Tin - (mt - m)is 



(30) 



(31) 



where sT = §T- 



and j £ = jef+jej.. Defining 
jheat = je — (5T/2)j„ [36], the Kubo formula gives us the 
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following spin Seebeck coefficient: 
1 



TSe = Pe 



dt[(J heat (0) ■ 3 sp in(t)) 



3VT J0 

- (ftf - /^)(J spm (0) • 3 spin (t))] (32) 
1 



WT 



dt(J heat (0) ■ 

u spin 

(*)> 



(33) 



where a s is the spin conductivity given in Eq. (2) and 
Mt ~~ = Tln(n^/n^). Therefore, in this representa- 
tion, the spin Seebeck coefficient (and thus also the spin 
Peltier coefficient) always carries an additional term of 
the spin conductivity a s multiplied by — ln(nf/nj_). This 
is the origin of that term in Pe' in Eqs. (15), (17), and 
(25) . We can understand the reason why the spin Seebeck 
coefficient Se' does not include such a term from the fol- 
lowing observation: At mechanical equilibrium and high 
temperature, we have 

V(Mf- W )=lnf5 t> ) VT + T^—Vp. (34) 



Therefore, that Vp = implies V(/Uf — /Ltj.) = 
ln(n-|-/n^)VT. The spin current now becomes 

jspin = -cr s V(^ t - ^) - 5eVT (35) 

= -D s nV P - (a, ln(n t /n 4 ) + Se) VT (36) 

= -D s nVp - Se'VT . (37) 

Thus the additional ~ \n{n^/n^) term from Eq. (34) 
exactly cancels the similar term from Eq. (33). Conse- 
quently, what we have computed for Se' in the previous 
section is the first term in Eq. (33) to which the sign 
argument in section IV should be applied. 

We can extract more information from Eq. (31). 
In jn = frame, the heat current is jheat — je — 
T\n{n^/n\)y S pin- Thus, heat conductivity and spin 
Peltier coefficient always have contributions coming from 
spin current. Therefore, in the format of Eq. (5) we can 
write 

k' = k[ - T\a{n t /n i )Se' (38) 
Pe' = Pe' l -T\n{n t /n i )D s , (39) 

where k\ and Pe\ are first terms which do not originate 
from spin current. Results in the previous section clearly 
show this. 

Lastly, we study the thermal conductivity, K, in Eq. 
(2). 

" dt<3 heat (t)-3 heat {0)> (40) 



WT 2 jq 
1 



WT 2 



dt < 3 heat (t) ■ 3 heat {0) > 



2Tln[ ^ \Se' 



T In 



(41) 



The above structure of thermal conductivity is explicitly 
shown in Appendix C. 



VIII. CONCLUSION AND OUTLOOK 



We have studied diffusive spin and heat transport in 
a two-species atomic Fermi gas with short-range inter- 
action at general polarization, temperature and scatter- 
ing length. Using "power-counting", we first estimated 
the spin and thermal diffusivities in all regimes. We 
suggested a method to measure the thermal diffusivity, 
which has not yet been experimentally measured. 

Our main focus was on the magnetocaloric effects, 
namely the spin Seebeck and spin Peltier effects. Ob- 
serving the connection between the interaction strength 
and the dependence of the scattering rate on the rela- 
tive momentum of two atoms, we were able to develop 
a qualitative argument for the signs of magnetocaloric 
effects. Near unitarity, magnetocaloric currents are in 
the same direction as the "primary" spin and heat cur- 
rents, while their directions are reversed as we move to far 
away from unitarity. We then quantitatively estimated 
diffusivities and magnetocaloric effects in the classical 
regime using approximate solutions of the Boltzmann ki- 
netic equation, thereby confirming the "power-counting" 
estimates of diffusivities and the sign argument for the 
magnetocaloric effects. We also proved the scaling of the 
transport coefficient is robust against to orders of approx- 
imation, in Appendix A. 

Remaining in the classical regime, we proposed an ex- 
perimental procedure to detect the spin Seebeck effect as 
a qualitative effect: a nonmonotonic relaxation towards 
equilibrium. This method is nice in that it does not re- 
quire knowledge of the equation of state of the gas. In 
order to confirm these results and obtain a better under- 
standing of the origin of these magnetocaloric effects, we 
also performed a controlled perturbation to the exactly 
solvable Maxwellian model. This approach agrees well 
with the approximate solutions of the Boltzmann equa- 
tion. 

To obtain the corresponding transport behavior quan- 
titatively in the quantum degenerate (but normal) 
regime, we need spin-polarized Fermi liquid theory. We 
hope to report on this soon in a forthcoming paper. 
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Appendix A: Approximate solution of Boltzmann 
equation to all orders 



Therefore, Eq. (13) simplifies to 



1. k' and Se' 

First let's consider k' and Se' when Vp = 0. 
Observe that Vp = and mechanical equilibrium con- 
dition imply 



Vji 

n 



(Al) 



4>i 



T 



(A2) 



where we abbreviated Cm + Cuj = C«- Following 
the procedure explained in the main text, we can ob- 
tain the following set of linearly independent equations 
(m = l,2,--- ,L). 



n T A 3 [A,C 0T + /3iC lt + • • • /? L C Lt ] + n ; A 3 [p C Qi + fcC n + ■■■ (3 L C Ll ] = 

"oo (Cot " Co;) + «oi (Ci-f - C14.) H h a L (Cl^ - C14) = A Q 

n^X 3 [a m0 (Cot ~ C ;) + (a m ifCit + OmU^u) + 1" (a m LtCLt + "mL^Ci;)] = A m 

n t A 3 [-a m0 (Cot ~~ Co;) + (aW4.Cit + a ml1 .C u ) H h (a mLi C L ^ + a mL ^C L ^)] = A m 



(A3) 
(A4) 
(A5) 
(A6) 



Here, {a 00 , a 0n , a m0 , a mna } (m,n = 1, 2, . . . , L and a =f 
, 4,) are dimensionlcss numbers which may depend on 
(A/a) 2 through the exponential integral and incomplete 
Gamma functions and implicitly depend on the order of 
approximation L but do not explicitly depend on temper- 
ature and densities. {At, fit} are determined by Gaussian 
integrals as following: 



A,= 



d 3 q 

(2tt) 3 

■)3+2£ A+< 



SI 

47T 

er 



(A7) 



d3q a 2e a 2 e- q2 / 47T = -2 3+2e ^ 
(2tt) 39 9z 3 



(A8) 



where T[a;] is the Gamma function and and are 
purely numerical numbers independent of physical pa- 
rameters. Note that the first equation, Eq. (A3) fixes 
the frame to be jt + jt = 0. From dimensional anal- 
ysis and the symmetry under t*^4-j we may write an 
ansatz solution of the above equations in the following 
form (m = 1,2, ... , L): 



Cot 
Cot 

a 



C m t 

where {at,bt,c Q } 



1 

A 3 " 
1 

A 3 " 

1 

A 3 " 

1 

A3 



nt 

a m 
n+ 

bin 

nt 
0,1,2, 



H 
ap_ 

H 

bm 
H 



, L) are 2L 



(A9) 
(A10) 
(All) 
(A12) 
3 un- 



known dimensionless numbers which may depend only 
on (A/a) 2 . Since the above set of equations should hold 
for any values of and nj,, once we insert the ansatz to 
the above equations, we obtain 2L + 3 linear equations 
in terms of {a^, 6^, Cq}. Since we are still using the same 
set of coefficients {aoo, aon, ctmo, ct mna -} and {fit}, it is 
easy to see that once the original set of equations is lin- 
early independent (which is a necessary condition to have 
an approximate solution of the Boltzmann equation) , the 
set of descendent equations is also linearly independent. 
Therefore, once we determine all {a^,6^,Co}, the above 
ansatz is the (L'th order) approximate solution of the 
Boltzmann equation. 

When Vp = 0, we know that j sp i n = —Se'VT. It is 
straightforward to show that 



Se' 




(A13) 



nt 



A) Co 



(A14) 



Since {ae,bt, Co} satisfies Eq. (A3), we have two identi- 
ties: 



1 

^2 (3 tat = --/3 c 

L 

Mi = ° • 



(A15) 
(A16) 



Inserting these identities to Eq. (A14), we finally obtain 
the full polarization and temperature dependence of the 
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Seebeck coefficient 



2?r). 



2. D s and Pe' 



Se = -c 



(A17) 



This proves that the scaling of Eq. (19) is indeed true 
at all orders of approximation. Furthermore, we see that 
the dimcnsionless scaled function h^(x) (L is the order 
of approximation) is 



h L ((\/a) 2 ) = ttcq 



(A18) 



It should be noted that cq implicitly depends on the order 
of approximation. 

For the thermal conductivity k', we are only interested 
in the first term, k'x, which directly comes from the energy 
current j e . 

1 h L 

K 'i = - ^ TO E (MntCtt + ni C n )) , (A19) 
n m t=o 

where Be = 47r(5/2 + t)p$_. Plugging in ansatz solution 
and using Eq. (A3), we obtain k^. 



£=1 



in 



A 3 



ntH 



Ifa (A20) 



Again, we emphasize that a a and bg are function of (A/a) 2 
and implicitly depend on the order of approximation. 
This proves that the algebraic form of the polarization 
dependence of the thermal conductivity obtained in the 
main text remains to all orders of approximation. 



Now we study D s and Pe' when VT = (thus Vn = 
0). As we did in the previous subsection, we want to 
express Vn| and Vrij in terms of riVp. 



m 



j riVp 
2m 



(A21) 



where = 1 and = —1. Then, the ansatz, Eq. (13), 
takes the following form (i ^ j): 



4>i 



2^1 n, n 



qf q, • Vp (A22) 
(A23) 



Once we define 



C/_i = 



n ( Cm C, 



n ; 



(A24) 



we reduce the system similar to the previous case. 
The 2{L + 1) linearly independent equations are (m = 
1,2,---L) 



n T A 3 l3 Q C ot + (3iC lt + ...(3 L C Lt - n ; A 3 f3 Q C ai + PiC n + . . . p L C Li 



"oo (pot + CoiJ + «oi (pit + CiiJ 



Cn + Cu) H h a 0L Cj* + C, 



%%A 3 
n 

i t n;A 3 



(pLt + 

ot m o (pot + Coij + ( atmit&it — a miiC\^j + • • ■ + (ot m i,^Ci,t — a m h\.Ch^ 

Ot m (pot + Co^j + (a m ilCi^ — a m i^Ci±j + • • • + (a m L^CLt ~ &mLtChi 






(A25) 


^0 


(A26) 


Am 


(A27) 


A 


(A28) 



a's and /3e are same as in the previous subsection and 
At = ft. 

Although coefficients of Eqs. (A27) and (A28) are lin- 
early dependent, once we combine them, we obtain an- 
other L linearly independent equations. 

L 

J2(a mn t - owX^W - ^W) = ■ ( A29 ) 

n=l 

Together with Eq. (A29), we have 2(L + 1) linearly in- 
dependent equations that uniquely determine all Cu. 



First observe that the solution of Eq. (A29) is trivial; 
C m t = C m i. Then, we use the following ansatz [37]: 



C, 



at 



Cot ~ 



C, 



mt 



1 


/ ao 




A 3 


\ n t 


H J 
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Ao 


ao ' 


A 3 


\ n t 


n i ) 


1 


TlC m 


— C m 


A 3 


ntni 



(A30) 

(A31) 
(A32) 
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{do, bo,c m } are dimensionless numbers that may only de- 
pend on (A/a) 2 . Substituting the above into original lin- 
ear equations, we obtain the same linearly independent 
equations in terms of {do, bo, c m }. Therefore, once we de- 
termine them, we have the approximate solution of the 
Boltzmann equation. 

Following the same procedure when we obtained Se' 
and k' 1; we can show that 

D, = — ^(So-So) (A33) 

Pp i = - 4s E - — - ■ ( A34 ) 

where Pe' x is the first term in the Peltier coefficient. Once 
we scale the scattering length with A, the polarization 
and temperature dependence of D a and Pe^ at an arbi- 
trary L remains the same as for L = 1. This proves that 
the scaling structure of the transport coefficients remains 
the same to all orders of approximation. 

Appendix B: Calculation of the spin Seebeck effect 

Here we present derivations of signatures of the spin 
Seebeck effect in detail. First, apply a long wavelength 
temperature modulation to the system. The mechanical 
equilibrium condition implies that a temperature modu- 
lation is equivalent to a total density modulation. Ini- 



tially the gas is uniformly polarized so that the density 
of each spin component is 

rii(t = 0, z) = riio + 5riio cos wz (Bl) 

Sn i0 = —5n , (B2) 
n 

where Sno is the maximum value of total density devia- 
tion from the average density. In order to calculate the 
change of densities of each species, which is directly mea- 
surable in experiment and contains the signature of the 
spin Seebeck effect, we need the space-time dependence of 
the density of each species. We will use the diffusion ma- 
trix and the continuity equation to derive the space-time 
dependence of densities. Since the continuity equation 
for the particle number density is dtfii + V ■ jj = 0, it 
is practical to write diffusion equation in the format of 
Eq. (1). With this initial condition, a nonuniform par- 
ticle current flows, but in this classical limit, there is no 
net energy current. Thus we work in the reference frame 
where the local energy current vanishes (j e ^ + j e i = 0). 
Thanks to Galilean invariance, this choice of a reference 
frame does not affect the physics. 

Following a similar procedure as described in the main 
text, but in this zero energy current frame, we can deter- 
mine all coefficients Cuj {t = 0, 1) and express the par- 
ticle current of each species in terms of Vn| and Vnj. 
Then, we can write currents in the diffusion matrix for- 
mat as in Eq. (1): 



Near unitarity, 



Far from unitarity, 



9tt 3 / 2 (h (T_ 
9728^2 1 m \T F 



3/2 N 



365n^+354n^n^-f 



733n l 



381n^+42n T n 4 +405nJ 




896\/2 \ m(k F a) 



nni 




405n^+42n T n4,+381nJ 


733n^+354u T n i +365nJ 






' 145n^ + 158n t n 4 +289nJ 


137n^-14n t ri4. + 125nJ 






125n^-14n t n 4 + 137raJ 


289n^ + 158n T ri4. + 145nJ 



nn-f 



nn-f 



/Vn t 
\^H 



\^H 



(B3) 



(B4) 



These matrix equations contain two eigenmodes, one is the thermal mode with eigenvalue St (j-j- and jj. are in the 
same directions) and the other is the spin mode with eigenvalue 23 s (j-f and jj, are in opposite directions). As expected, 
Di > S) s for all polarizations in both of these limits. Applying continuity equations to the above diffusion matrix 
equations while keeping all differential operators linear (we restrict ourselves in a linear response theory), we obtain 
two-component diffusion equations: 



0_ 
di 



D 



tt 



n 



v 2 



Therefore, we can immediately write the time evolution of each species in terms of the eigenmodes as 



n t (t,z) 

h(^ z ) 



n T0 
Ho 



8 Ho 



cos wz 



-f3e 



-S) T w 2 t 



$Ho 



(B5) 



(B6) 



(7, 1) and (£, 1) are the eigenvectors of the spin mode and the thermal mode, respectively. 7 is negative while £ is 
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positive, a and /3 are determined by the initial condition, Away from unitarity, 



(B7) 



It is straightforward to derive 7, £, a, /3, S) s and £>t- For 
example, 



1 / «ot > 

a = C 

7 - C 







1 



(B8) 
(B9) 



Explicit expressions for 7, £, S) s and £>t are fairly lengthy 
so they are omitted here. 

The spin Seebeck effect is most apparent when observ- 
ing the densities at the edges of the system (wz = or 
wz — 7r) so we choose to focus on the cold edge, z = 0. 
The uniform pressure condition implies that both major- 
ity and minority densities are initially higher than aver- 
age at the cold side. Thus, initially Sn^(t = 0, z = 0) > 0. 
Even though the initial condition can be chosen to be 
the same for both unitarity and far away from unitarity, 
the time evolution of density of each species is qualita- 
tively different for these two limiting cases. As argued 
in the main text, the minority density deviation has an 
extremum before it relaxes to zero at unitarity whereas 
the majority density deviation has an extremum far away 
from unitarity. Therefore, a change in the sign of a cer- 
tain species is a unique signature of the spin Seebeck ef- 
fect. We can express the magnitude of the spin Seebeck 
signal in dimensionless form by ^y^p- This quantity 
is plotted in figures 3 and 4. At time ti^ ext (i =-1 near 
unitarity and i —'f away from unitarity), the density of 
the species i reaches an extremum. It is easy to derive 
formal expressions of ti^ ex t and a dimensionless measure, 



Sn i0 

Near unitarity, 



1 



(® T -D s )w 2 



In 



= a exp 
+ /9exp 



-In 
In 



D T /3 



(BIO) 



(Bll) 



Sn^o 



+ 



°7 

(K 

«7 + /3C 



exp 
exp 



S s a7 



(B12) 



In 

ln 



S s ck7 
(B13) 



.0) 



Srii, 



is plotted in figure 6. As expected, the 

unitarity limit shows a stronger signal of the spin Seebeck 
effect than far away from unitarity. 

Another consequence of the spin Seebeck effect is the 
change of polarization. From Eqs. (B2) and (B6), we can 
express the polarization as a function of time. Keeping 
only linear term of Suq, we obtain 



6p(t,0)=p(t,0)-po 



Sno 
-Po 



Ho 
n 



("(7 
(a(7 4 



1 e 



+ /3(C-l)e 



De 



(B14) 



where po = (n^o — ^4.0 ) /«o • Figure 5 is the plot of Sp(t , 0) 
normalized by dno/n a . 



Appendix C: Manifestation of the Onsager relation 

For completeness, we present diffusion matrices of ap- 
proximate solutions and the solution of the first order 
perturbation of the Maxwellian model in the form of Eq. 
(2) where the Onsager relation should be explicit. We 
simply transform the set of driving forces (VT, nVp) to 
another set of driving forces (VT, V(^-j- — /ij,)) associated 
with the entropy production. 



For approximate solutions we have the following: 
At unitarity (|a| >• A), 



iheat 
\ispin , 



45tt 3 / 2 
608^ 



T 
Tf 



3/2N 



T 5nT 



In 



39n-j-n^ 
5nT 



VT 

v(/xt - m) 



(CI) 
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where k u is 



5 / (73rd + 82n tni + 73nf) 

n u — — n 

16 \ n t n; 



In 



n^-\ ( 10(n| - u|) - 39n t n;ln(n t /n;)^ 



5n 



Far away from unitarity (|a| <C A), 



J/ieai 
J sp in 



where K a is 



15tt 3 / 2 
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-(n t -n 4 )-i 3 ^±ln(2l) S 
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For the first order perturbation of the Maxwellian model, we have the following: 



iheat 
J spin 



where Km is 
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107T6 J {n^ 
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The above matrices clearly exhibit the Onsager rela- section VI. Lastly, the structure of thermal conductivity 
tion. Also each off-diagonal element carries a term pro- (Eq. (41)) is evident, 
portional to — a s ln(n+/ni) as mentioned in Eq. (33) in 
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